Weakly symmetric stress equilibration for hyperelastic materialmodels

A stress equilibration procedure for hyperelastic material models is proposed andanalyzed in this paper. Based on the displacement-pressure approximation computed with a stable finite element pair, it constructs, in a vertex-patch-wise manner, an $H(div)$-conforming approximation to the first Piola-Kirchhoff stress. This is done in such a way that its associated Cauchy stress is weakly symmetric in the sense that its anti-symmetric part is zero tested against continuous piecewise linear functions. Our main result is the identification of the subspace of test functions perpendicular to the range of the local equilibration system on each patch which turn out to be rigid body modes associated with the current configuration. Momentum balance properties are investigated analytically and numerically and the resulting stress reconstruction is shown to provide improved results for surface traction forces by computational experiments

Weakly symmetric stress equilibration and a posteriori error estimation for linear elasticity

A stress equilibration procedure for linear elasticity is proposed and analyzed in this paper with emphasis on the behavior for (nearly) incompressible materials. Based on the displacement-pressure approximation computed with a stable finite element pair, it constructs an $H (\text{div})$-conforming, weakly symmetric stress reconstruction. Our focus is on the Taylor-Hood combination of continuous finite element spaces of polynomial degrees $k+1$ and $k$ for the displacement and the pressure, respectively. Our construction leads then to reconstructed stresses by Raviart-Thomas elements of degree $k$ which are weakly symmetric in the sense that its anti-symmetric part is zero tested against continuous piecewise polynomial functions of degree $k$. The computation is performed locally on a set of vertex patches covering the computational domain in the spirit of equilibration \cite{BraSch:08}. Due to the weak symmetry constraint, the local problems need to satisfy consistency conditions associated with all rigid body modes, in contrast to the case of Poisson's equation where only the constant modes are involved. The resulting error estimator is shown to constitute a guaranteed upper bound for the error with a constant that depends only on the shape regularity of the triangulation. Local efficiency, uniformly in the incompressible limit, is deduced from the upper bound by the residual error estimator

Weakly symmetric stress equilibration and a posteriori error estimation for linear elasticity

This paper proposes and analyzes a posteriori error estimator based on stress equilibration for linear elasticity with emphasis on the behavior for (nearly) incompressible materials. It is based on an H(div)-conforming, weakly symmetric stress reconstruction from the displacement-pressure approximation computed with a stable finite element pair. Our focus is on the Taylor-Hood combination of continuous finite element spaces of polynomial degrees k + 1 and k for the displacement and the pressure, respectively. This weak symmetry allows us to prove that the resulting error estimator constitutes a guaranteed upper bound for the error with a constant that depends only on local constants associated with the patches and thus on the shape regularity of the triangulation. It does not involve global constants like those from Korn's in equality which may become very large depending on the location and type of the boundary conditions. Local efficiency, also uniformly in the incompressible limit, is deduced from the upper bound by the residual error estimator. Numerical results for the popular Cook's membrane test problem confirm the theoretical predictions

Sea Ice conditions within the Antarctic Marginal Ice Zone in winter 2017, onboard the SA Agulhas II

Our knowledge of sea ice variability, which contributes to the detection of the Antarctic climate change trends, stems primarily from remotely sensed information. However, sea ice in the Southern Ocean is characterized by large variability that remains unresolved and limits our confidence on the remotely sensed products. Therefore, the in situ sea ice observations presented (according to the ASPeCt protocol) provide a greater understanding of the Antarctic sea ice environment - on a local scale - and allows us to evaluate remotely sensed products